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Abstract 

In the framework of f{T) theories of gravity, we solve the field equations for /(T) = T in the 
weak-field approximation and for spherical symmetry spacetime. Since f(T) = T corresponds to 
Teleparallel Gravity, which is equivalent to General Relativity, the non linearity of the Lagrangian 
are expected to produce perturbations of the general relativistic solutions, parameterized by a. 
Hence, we use the f{T) solutions to model the gravitational field of the Sun and exploit data from 
accurate radio-tracking of spacecrafts orbiting Mercury and Saturn to infer preliminary bounds on 
the model parameter a and on the cosmological constant A. 
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I. INTRODUCTION 


The geometrical structure of Teleparallel Gravity (TEGR) [l-3| is constituted by a 
Riemann-Gartan space-time, endowed with the Weitzenbock connection; the latter, unlike 
the Levi-Givita’s one, is not symmetric and gives rise to torsion but it is curvature-free. On 
the other hand, the space-time of General Relativity (GR) is endowed with the symmetric 
Levi-Givita connection, which is torsion-free and gives rise to the Riemannian curvature. 
Space-time curvature defines the gravitational interaction in GR: the field equations are 
obtained from the Einstein-Hilbert Lagrangian, which contains the curvature scalar R. By 
contrast, in TEGR torsion plays the role of curvature, and the field equations are obtained 
from a Lagrangian containing the torsion scalar T. Even though the two theories have a dif¬ 
ferent geometrical structure, they share the same dynamics: this means that every solution 
of GR is also solution of TEGR. 

The interest in gravity theories with torsion dates back to Einstein: in fact, in his attempt 
to formulate a unihed theory of gravitation and electromagnetism, Einstein made use of 
tetrads and torsionjd-Gl. There is another important approach in which torsion is used 
to determine the geometry of space-time, the so-called Einstein-Gartan Theory (EGT) of 
gravity: from a geometrical viewpoint, EGT is an extension of GR to Riemann-Gartan 
spaces, which have both curvature and torsion; in particular, torsion is related to the density 
of spin j?-^. Actually, the effects of spin and torsion are signihcant only at very high 
densities of matter: as a consequence, EGT is a viable theory of gravity since its observational 
predictions are in perfect agreement with the classical tests of GR; moreover, it is suggested 


that EGT will prove to be a better classical limit of a future theory of quantum gravity 


10 |. 


In Teleparallel Gravity the tetrads field plays the role of the dynamical field instead of 


the metric field. Recently [h|, the properties of Weitzen 


□oc 


with reference to the nonlocal generalization of GR 


12 


^ connection have been studied 


are a generalization of TEGR, in which the gravitational Lagrangian is an a na 
of the torsion scalar. Actually, these theories are not equivalent to GR 


15| . The so-called /(T) theories 
ytic function 
, 17|: this is 


why they 


puzzle 


18 


rave been considered in cosmology as candidates to solve the cosmic acceleration 


29|. Beyond cosmology, other aspects of f(T) theories have been investigated. 


such as exact solutions, in 


electromagnetic field 30l-l42 1. 


ower or higher dimensional space-time, also in presence of the 
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It is important to emphasize that f{T) have additional degrees of freedom 43|, [4^ with 
respect to TEGR (and, hence, to GR), that are exploited in the aforementioned cosmolog¬ 
ical applications; however, these degrees of freedom do not contribute to t he g ravitational 


radiation when considering the hrst order perturbation theory, as shown in 


45|. 


They are related to the fact that the equations of motion are not invariant under local 
Lorentz transformations: as a consequence, there exists a preferential global reference frame 
dehned by the autoparallel curves of the manifold that solve the equations of motion. So, 
even though the form of the space-time metric can be obtained on the basis of the physical 
symmetries, it is necessary to pay attention to correctly dehne the form of tetrad field. In 
particular, as clearly discussed in j^, a diagonal tetrad is not the best choice for non-flat ho¬ 
mogenous and isotropic cosmologies (Friedman-Lemaitre-Robertson-Walker universes) and 
for spherically symmetric space-times (Schwarzschild or Schwarzschild-de Sitter solutions). 
Spherically symmetric solutions in f(T) gravity are important also because these solutions, 
describing the gravitational held of point-like sources, can be used to constrain these theories 
with planetary motions in the Solar System. A weak-held solution, suitable to model the 
gravitational interaction in the Solar System, was obtained in for a Lagrangian in the 
form /(T) = T + aT'^, using a diagonal tetrad; a is a small constant which parameterizes 
the departure from GR. Observations of the rates of change of perihelia allowed to constrain 


the theory parameter a; tighter constraints on a were obtained in 


48| . It is interesting to 


point out that such a form of Lagrangian can be found also in Born-Infeld gravity 491], and 
its cosmological implications have been investigated in 50]. 

In a recent paper Mlj we studied spherically symmetric solutions, in the weak-held ap¬ 
proximation, for Lagrangians in the form /(T) = T + aT'^, with \n\ ^ 1, and solved the held 
equations using a non diagonal tetrad. We showed that, to lowest approximation order, the 
perturbations of the corresponding GR solutions (Schwarzchild or Schwarzchild- de Sitter) 
are in the form of power laws oc a and we estimated the a parameter for n = 2. 

In this paper, we use the solutions found in ^ to obtain new bounds on a and on the 
cosmological constant A, by exploiting data from tracking of spacecraft orbiting some of the 
major bodies of the Solar System like MESSENGER (Mercury) and Gassini (Saturn). 

This work is organized as follows: in SectionlTTlwe briehy review the theoretical framework 
of /(T) gravity and the held equations and, then, we give the solutions for spherically 
symmetric space-times, in weak-held approximation; in Section IIIII we infer bound on the 
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model parameters; eventually, summary and conclusions are Section IIVI 

II. WEAK-FIELD SPHERICALLY SYMMETRIC SOLUTIONS FOR /(T) GRAV¬ 
ITY 

In /(T) gravity the metric tensor can be viewed as a subsidiary field, and the tetrad field 
is the dynamical object. Given a coordinate basis, the components of the tetrad field e“ are 
related to the metric tensor by 

9i,y{x) = 7]ahel{x)el{x) , ( 1 ) 

where rjah = diag(l, —1, —1, —1). Notice that latin indexes refer to the tangent space while 
greek indexes label coordinates on the manifold. The field equations are obtained from the 
action^ 


S=j^Jf(T)ed*x + SM, ( 2 ) 

where e = det e“ = a/— det(^^ and Sm is the action for the matter fields. In Eq. ([2]) /(T) 
is a differentiable function of the torsion scalar T, which is defined as 

T = ’ (3) 

where the contorsion tensor is defined by 

SP =-(TP -T P + T P^ + -5PT ^ - -5PT ^ f4j 

^ ^ \ fii/ ' ufi J ' ^ av afi ’ V / 

and the torsion tensor is 

rp = - 9X) ■ (8) 

The variation of the action with respect to e“(x) gives the following field equations: 

e-'d,(e e,>S;-)fT + S + eJ‘S/''d,(T)fjr + \e'^f = ( 6 ) 

where 7^^ is the matter energy-momentum tensor and subscripts T denote differentiation 
with respect to T. 

^ We use units such as c = 1. 
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We are interested in spherically symmetric solutions of the held equations, so we write 
the space-time metric in the general form 


ds^ = - r^d^f , 


(7) 


where dVL^ = dO'^ + ?Av?6d(f)^. We follow the approach described in ^ to derive the held 
equations for the non diagonal tetrad: 


e 


a 


/ 0 0 ^ \ 
0 sin 9 cos 0 sin 6 sin 0 cos 9 

0 —r cos 9 cos (j) —r cos 9 sin cj) r sin 9 

y 0 r sin 9 sin 0 —r sin 9 cos cj) 0 J 


( 8 ) 


We do emphasize that, because of the lack of local Lorentz invariance, tetrads connected by 
local Lorentz transformations lead to the same metric - i.e. the same causal structure - but 
diherent equations of motions, thus physically inequivalent solutions. In particular, even 
in an extremely symmetric case like the static spherically symmetric one, a diagonal tetrad 
that gives back the metric in eq. ([7]) is not a good choice since the equations of motion for 
such a choice would constrain a priori the form of the Lagrangian. Actually, this is not the 
case for the non-diagonal tetrad chosen in eq. (IH]). In this case, the torsion scalar turns out 
to be 


2p-^h) (\ 4- eB{r)/2\ 

T = -[1 + ^ rA\r)] . 


(9) 


Exact solutions in vacuum and in presence of a cosmological constant of the above held 
equations (jb]) are discussed in in details. If a vacuum solution of f(T) has T = 0 
it will be a solution of TEGR (and hence GR) as well j^, 38|. In fact, the equations 
of motion with T = 0 reduce to the Einstein equations of motion with a cosmological 
constant A = and one can even avoid the Lambda term by choosing /(0)=0 and 

f'{0) 7^ 0. Here, we are interested in weak-held solutions with non constant torsion scalar, 
i.e. T' = dT/dr ^ 0: as shown in in this case the most general vacuum solution is 
diherent from the Scwharzschild - de Sitter one. This is interesting because new features 
arise with respect to the vacuum solutions of GR, so deviations from GR can be in principle 
observed. 

In the weak-held limit, that is suitable to describe the gravitational held in the solar 
system, we may write = 1-|-A(r), = 1-|-i?(r); furthermore, we conhne ourselves 
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to linear perturbations. As shown in Ml], if we consider Lagrangians in the form /(T) = 
T + where a is a small constant, parameterizing the departure of these theories from 
GR, and |n| 7 ^ 1, in presence of a cosmological constant A, we obtain the following solutions 


, 2GM r 
A{r) = - a 


2—2rx 1 

-- -Ar^ 

2n — 3 3 


( 10 ) 


2GM .^2-2n 1 

B{r) =-h OL- -(- 377 , + 1 + 2r?) + -Ar^ (11) 

V ZiTl O O 

for the gravitational held around a point-like source M. In particular, if A = 0, we obtain 
vacuum solutions. We see that the weak-held limit of Schwarzschild - de Sitter solution is 
perturbed by terms that are proportional to a. Moreover, the torsion scalar is 


T(r) = ^ + 2ar-^^2^^ {n + 1) 


( 12 ) 


For n = 2, the space-time metric is 

2GM 


ds^= 1 


-32^]dR- {1 + 


2GM ® \ 7 2 2 7r^2 

-h 96— I dr^ — r^dVL 


(13) 


Notice that the perturbations due to the non linearity of the Lagrangian have the same 
functional behaviour as in 47||, but different numerical coefficients: this is not surprising, 
since the held equations are not the same. 

On writing the perturbation terms due to the deviation from GR in the form 


Aa{r) = aanr'^ 5„(r) = abnv'^ (14) 

2^Sn—l 2^3n—l 

where a„ = -, = - (2n^ — 3n -|- l), we see that they go to zero both when 

2n — 3 2n — 3 

r —)■ cxo with n > 1 and when r —)■ 0 with n < 1. In the latter case, in order the keep the 
perturbative approach self-consistent, a maximum value of r must be dehned to consider 
these terms as perturbations of the hat space-time background. 

Eventually, we emphasize that our approach can be applied to an arbitrary polynomial 
correction p{T) to the torsion scalar, i.e. to /(T) = T + p(T): in particular, by writing 
an arbitrary function as a suitable power series, it is possible to evaluate its impact as a 
perturbation of the weak-held spherically symmetric solution in GR: the n-th term of the 
series gives a contribution proportional to 
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III. PHENOMENOLOGICAL BOUNDS ON a AND A FROM SOLAR SYSTEM’S 
ORBITAL MOTIONS 


As we have seen above, for Lagrangians in the form /(T) = T + aT” the perturbations 
to the GR solution are in the form ~ We focus on n > 1: actually, since realistic 

deviations from GR are expected to be small in the Solar System, we can consider, to a 
hrst approximation, the solution for n = 2. Moreover, we do include the contribution due 
the cosmological constant A: in other words, we consider a space-time metric in the form 


(IT^ . Hence, here we will phenomenologically infer preliminary bonne 


latest results in the held of the Solar System planetary ephemerides 


52 


s on a, A from the 


55l |. In doing so, we 


will suppose that orbital motions follow metric geodesics. Indeed, it is useful to remember 
that in theories with torsion there is a distinction between the test particles trajectories: 
autoparallels, or affine geodesics, are curves along which the velocity vector is transported 
parallel to itself, by the space-time connection; extremals, or metric geodesics, are curves 
of extremal space-time interval with respect to the space-time metric 56(]. While in GR 
autoparallels and extremals curves do coincide and we can simply speak of geodesics, the 
same is not true when torsion is present. 

To infer bounds on a, A, we will look at the geocentric right ascension RA, declina¬ 
tion DEG and range p, which are directly observable quantities usually measured in ac¬ 
curate tracking of spacecraft orbiting some of the major bodies of the Solar System like 
MESSENGER (Mercury) and Gassini (Saturn). We will exploit the existing time series 
ARAexp(t), ADEGexp(^), Apexp(^) of the post-ht residuals of such observables produced by 
dedicated teams of astronomers by comparing long data records with the corresponding the¬ 
oretically computed observations on the basis of accurate multiparameteric models of the 
orbital dynamics of the probes themselves, of the planets, of the propagation of the elec¬ 
tromagnetic waves, of the measuring devices, etc. In such data reductions, dozens of such 
parameters like, e.g. the planetary initial conditions, etc., are determined in a least-square 
sense. 

Usually, only standard dynamics in terms of the currently accepted Einsteinian laws of 
gravitation is modeled, so that, in principle, the residuals ARAexp(t), ADEGexp(t), Apexp(^), 
which, by construction, account for any unmodeled and mismodeled effects, should be 
affected also by the modihed gravity studied here. Thus, we will produce time series 
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ARA(t), ADEC(t), Ap(t) of simulated residuals due to given values of a and/or A and com¬ 
pare them to the real residuals existing in the literature to infer the limiting values of a, A 
which makes our theoretical signatures compatible with ARAexp(t), ADECexp(t), Apexp(^)- 
Although not strictly rigorous, our method is able to yield reasonable answers about what 
can be inferred on the basis of the state-of-the-art of the planetary ephemerides. Indeed, 
if, on the one hand, the exotic dynamical models one is interested in should be explicitly 
modeled and dedicated solved-for parameters should be determined in ad-hoc newly gen¬ 
erated ephemerides re-processing the same data records, on the other hand the fact that 
our model is relatively simple, implying just one/two extra-parameter(s), helps in assuring 
that the results we obtain are plausible. Moreover, it is just the case to remark that our 
approach, applied by a plethora of other researchers to different observation-related quan¬ 
tities determined for different scopes with respect to the subsequently arisen goals of test¬ 
ing/constraining this or that alternative modihed gravity, has always been widely followed 
in the literature so far in a number of scenarios, including also, e.g., the PPN parameters, 
the Pioneer Anomaly, etc. In some cases, like MOND and the Pioneer Anomaly 

58l-l6l|. both approaches were applied giving substantially equivalent results. 

In order to produce our simulated residuals, we hrst numerically integrate the Solar 
System’s barycentric equations of motion of the Sun, the eight planets and the dwarf planet 
Pluto with and without the a/A additional accelerations. Both the integrations share the 
same initial conditions, retrieved from the HORIZONS web interface maintained by NASA 
JPL. Then, for a given Solar System’s major body and for suitably chosen values of a 
and A, we take the difference between both integrations and numerically calculate the time 
series of its RA, DEC, range which show the characteristic pattern of the impact of the 
modihed gravity considered on such observables: they are our theoretical templates, to be 
compared with the real residuals. In view of the different dependence on r of the a and A 
extra-accelerations, we will use Saturn for A (Figure [1]) and Mercury for a (Figure [2]). 

The orbit of Saturn was recently constrained to a very accurate level from the analysis 
of long records of various types of observations from the Cassini spacecraft, engaged in the 
exploration of the Kronian system since 2004: its end is currently scheduled for^ late 2017, 
with the atmospheric entry into the planet. The geocentric range of the ringed planet is 


^ See http://saturn.jpl.nasa.gov/mission/satumtourdates/saturntourdates2017/ on the WEB. 










currently known to a 0.1 km level ^ from the precise radio-tracking of the spacecraft, 
while its RA and DEC residuals were recently reduced down to a few milliarcseconds level 
0 by means of the Very Long Baseline Array (VLBA) technique. In Figure [H we show 
our simulated signatures for 

A = 5 X 10"^^ m-2. (15) 

It turns out that the range is more effective in constraining A since a larger value of it would 
compromise the compatibility of our Ap{t) signature with the Cassini-based range residuals. 
The resulting bound on A of f[T5D is one-two orders of magnitude tighter than those inferred, 


within different theoretical and phenomenological frameworks, in, e.g.. 


M 




^1 


As far as Mercury is concerned, its trajectory is nowadays mainly determined by range 
measurements to the MESSENGER spacecraft after its orbit insertion in early 2011; the 
mission ended this year, with a guided crash of the probe into the planet’s surface. Figure 
5 of depicts the range residuals covering about 6 months produced at JPL, NASA 
with a recent version of the DE ephemerides; they appear to be far smaller than 50 m. In 


53l |. an independent team from IMCCE processed one and half years of radioscience data 


of MESSENGER with the INPOPlSa ephemerides by constraining the geocentric range 
residuals of Mercury down to |Ap|g^p < 8.4 m. In Figure [2l our «-induced simulated range 
residuals obtained for 

a = 5 X 10"^ m^ (16) 

are shown; a larger value for a would make them incompatible with the existing 
MESSENEGR-based real residuals in I53II . The bound of flTBD is about one-three orders 


of magnitude better than those in 


5l| . As remarked, e.g., in |^, actua 


be even tighter if a smaller value for | Ap|exp was to be used. Indeed, neither in 52| nor in 53| 


ly, it could 


the Solar gravitomagnetic held was modeled, so that it is likely that most of the Hermean 
residuals produced so far are due to such general relativistic dynamical effect, which should 
be subtracted by leaving an even smaller unmodelled/mismodeled signature. However, we 
conservatively assume (ITHll . 

IV. SUMMARY AND CONCLUSIONS 


In this paper, we used the latest results from Solar System planetary ephemerides to study 
the impact of f{T) gravity. In particular, we focused on spherically symmetric solutions 
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Saturn: right ascension 



Saturn: declination 




2004 2006 2008 2010 2012 2014 


t (yr) 


FIG. 1: Numerically produced time series of RA, DEC, p of Saturn perturbed by A = 5 x 10 
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m“^. For each of the three Kronian observables considered, they were calculated as differences 
between two numerical integrations of the SSB barycentric equations of motion of the Sun, its 
eight planets and the dwarf planet Pluto from 1914 to 2014 with and without the A term. Both 
integrations shared the same initial conditions in rectangular Cartesian coordinates, retrieved from 
the HORIZONS web interface, and the same standard dynamical models, apart from A itself. 
Thus, such curves represent the expected A-induced signatures ARA, ADEC, A^. The patterns 
and the size of the present signals can be compared with the range residuals by 
km) and those for RA, DEC by 54] (|ARA|g^p , jADECIg^^ < 4 milliarcseconds). 


m (|Ap|e.p < 0.1 


I exp 
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Mercury: range 



FIG. 2; Numerically produced time series of p of Mercury perturbed by a = 5 x 10“^ m^. It was 
calculated as the difference between two numerical integrations of the SSB barycentric equations 
of motion of the Sun, its eight planets and the dwarf planet Pluto from 1914 to 2014 with and 
without the a term. Both integrations shared the same initial conditions in rectangular Carte¬ 
sian coordinates, retrieved from the HORIZONS web interface, and the same standard dynamical 
models, apart from a itself. Thus, such curve represents the expected a-induced signature Ap. 


The pattern and the size of the present signal can be compared with the range residuals by 


^1 


(lApIg^p < 8.4 m). Note that most of such a figure may be due to the unmodeled Lense-Thirring 
effect, which should be subtracted by leaving a smaller amplitude for the Mercury’s reasiduals. 


obtained for Lagrangians in the form /(T) = T + aT^. Since the case f{T) = T corresponds 
to TEGR gravity, which is equivalent to GR, a is a constant that parameterizes the departure 
of these theories from GR, hence it is expected to be small. Due to the lack of Lorentz 
invariance for local transformations, we used a non diagonal tetrad to obtain the solutions 
for the held equations and we also considered the presence of a cosmological constant A. 

The presence of the cosmological constant naturally suggests the comparison of the con¬ 
straints obtained in this context with those coming from cosmological solutions. Indeed, an 
alternative model of gravity should be phenomenologically viable and congruent in different 
settings and at least reproduce the amazing coherence of General Relativity. For the f{T) 
model a complete study of exact cosmological solutions coming from power law model via 
the Noether symmetry approach is presented in 63|], while the comparison with cosmological 
data has been reported in ISj] and the imprints on galaxy clustering and weak gravitational 
lensing are derived in 29|]. Nevertheless, a direct confrontation of the results is not possible, 
as it can be deduced by a rapid inspection of the torsion scalar, that is negative dehnite for 
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the diagonal tetrads suitable for the solutions in the cosmological setting ( with no spatial 
curvature). The correct way to proceed would be to derive the junction conditions for the 
solutions describing the different regions of spacetime one wants to match, to make sure of 
their accomplishment and then to infer something on the parameters. At the present stage 
the junction conditions for the f{T) model have been derived in 6^, and they turn out to 


be more restrictive than in GR since they depend both on the detail of the model considered 
and on the choice of the tetrads. Actually, it plays a crucial role due to the fact that the 
induced tetrads on the boundary have to be continuous; we are going to investigate this 
issue in future works. 

We used the aforementioned spherically symmetric solutions to describe the Sun’s grav¬ 
itational held and used the Solar System data to constrain the model parameters a and 
A, for the case n = 2, since for n > 2 the impact of the non linearity of the Lagrangian 
is expected to be smaller. To this end, we considered the right ascension, declination and 
range measured for the spacecrafts MESSENGER and Cassini, orbiting Mercury and Sat¬ 
urn, respectively. Our approach is based on the employment of the time series of the post-ht 
residuals of these observables: in fact, since in the models of orbital dynamics only standard 
GR is taken into account, deviations from GR should affect the residuals. So, we calculated 
time series of simulated residuals due to given values of a and A, and then we compared 
them to the actual values of the residuals to obtain the limiting values of these parameters 
that are compatible with the observations. This kind of approach, even though is not strictly 
rigorous, can be anyway used to infer plausible results, as we discussed in Section IIIIl 

As a result, the bounds that we obtained both for the cosmological constant A and for 
the parameter a, are tighter than those available in literature deriving by Solar System 


dynamics. In particular, for the cosmological constant the bound is A = 5 x 10 
while for the /(T) Lagrangian parameter, the bound is a = 5 x 10“^ m^. 


-44 j^-2 
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